Asymptotic shape of the region visited by an Eulerian Walker 
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We study an Eulerian walker on a square lattice, starting from an initial randomly oriented 
background using Monte Carlo simulations. We present evidence that, for large number of steps N, 
the asymptotic shape of the set of sites visited by the walker is a perfect circle. The radius of the 
circle increases as TV 1 / 3 , for large TV, and the width of the boundary region grows as N a/3 , with 
a = 0.40 ± .06. If we introduce stochasticity in the evolution rules, the mean square displacement 
of the walker, (R%) ~ TV 2 ", shows a crossover from the Eulerian (y = 1/3) to a simple random walk 
{y — 1/2) behaviour. 

PACS numbers: 05.50.+q, 02.50.-r 



I. INTRODUCTION 

There has been a lot of interest in the study of asymp- 
totic shape of growing clusters. One of the earliest of such 
studies was by Richardson, who showed that the asymp- 
totic shape of the infected region in an epidemic model 
has linear segments [l|. In the Eden model Q, which 
models an epidemic without recovery, it has been shown 
that the asymptotic shape of the growing cluster is not a 
perfect circle Q . The shape of growing clusters has also 
been studied in sandpile models. In the abelian sandpile 
model, in two dimensions in an initial background of h 
particles at each lattice site [3, [f| , it was found that the 
cluster in general has a convex asymptotic shape, which 
becomes more circular as h is decreased, and tends to a 
perfect circle as h — > — oo. 

In this paper, we study the Eulerian walker (EW) 
model on the square lattice. This model is related to the 
sandpile model and was initially introduced by Priczzhcv 
et al. [1, 01 as a simpler variant of the sandpile model of 
self-organized criticality (SOC) Q. It has subsequently 
found applications in design of derandomized simulations 
of Markov chains Q, efficient information transfer pro- 
tocols in computer networks (lfj| . and modelling coevo- 
lution of virus and immune systems [ll[ . We study the 
model, starting from a disordered background, by Monte 
Carlo simulations. We have a single walker that moves on 
the lattice and we look at the shape of region visited by it 
which grows with the length N of the walk. Interestingly, 
we find evidence that this region is asymptotically a per- 
fect circle. The circular shape is reminiscent of the circu- 
lar shape in the rotor-router aggregation model studied 
by Propp (see [l^, Ell), where, for a special initial con- 
figuration, the region is almost a perfect circle [l4| with 
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departures from the circle being of order 1. The circular 
shape for the EW cluster is not very evident for small 
walks. For example, in Fig. [1] we have shown clusters 
formed by the EW of N = 10 5 and 10 7 steps. Clearly, 
only for large N, does the circular shape start to emerge, 
and it requires careful statistical analysis to see this when 
TV is not so large. 

The EW model can also be looked upon as a particular 




FIG. 1: (Color online) Cluster formed by an EW after 10 5 
(top) and 10 7 (bottom) steps. A site is assigned a color de- 
pending on whether the numbers of visits to it is in the range 
1-4,5-8 etc. for 10 s , and 1 - 16, 17 - 32 etc. for 10 7 
steps (key panel displayed on the right). The diameters of 
the clusters are 134 and 600 units respectively. 
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FIG. 2: (Color online) A portion of the square lattice. At 
each lattice site, the outgoing bond is shown by the arrow. 



limit of a growing self repelling walk, in which the walker 
preferably jumps along a bond which has been visited 
least number of times so far. This model was studied in 
Id by Toth and Veto In the zero temperature limit, 
in one specific variant, this becomes the EW model and 
a finite temperature corresponds to noise. It was found 
that the number of visits at a distance y from the origin 
satisfies the scaling function F(y) = 1 — y, for < y < 
1. We find that the same scaling function holds even in 
2d. We also study the model in the presence of noise, 
where there is a small probability e that walker goes in a 
direction not given by the EW rule. We find that a small 
noise changes the asymptotic behaviour. The diameter 
of growing region scales as A(e)N 1 ^ 2 in the presence of 
noise and as TV 1 / 3 in its absence. 

The paper is organized as follows. In Sec. [ill we define 
our model. In Sec. IIII1 we give details of the simulation 
without noise (i.e. e = 0). We find that the asymptotic 
shape of lines of constant average number of visits by 
an EW are perfectly circular, within statistical errors. 
We will argue that the variance of average number of 
visit at any distance from the origin tends to a finite 
number for large N and also obtain the scaling function 
for the average number of visits. In Sec. IIVI we discuss 
the case with noise, and finally we summarize our results 
in Sec. El 



II. MODEL 

The Eulerian walker is defined as follows: We consider 
a square lattice. We associate with each site an arrow 
which can point to along one of the four directions, de- 
noted by N, E, S and W (Fig. [5]). In the initial configura- 
tion, the direction of the arrow at each site is chosen inde- 
pendently, and with equal probability. We put a walker 
at the origin which moves on the lattice. The motion of 
the walker is affected by configuration of arrows on the 
lattice, which in turn affects the arrow configuration on 
the lattice. The walker follows the following rule: at each 
time step, the walker after arriving at a site rotates the 
arrow at that site in a clockwise direction by 90°, and 
then moves one step along the new arrow direction. 

It was shown in Ref. Q that on any finite graph, us- 




FIG. 3: (Color online) Mean square displacement B? N of the 
walker from the origin as a function of N for various e values. 
For the pure case (i.e. e = 0) -R^ grows as N 2 " with v — 1/3. 
In the presence of noise, there is a cross over from the Eulerian 
[y = 1/3) to a simple random walk behaviour (y = 1/2). 



ing the above rules, the walker eventually visits all sites 
and settles into a limit cycle which is an Eulerian circuit 
visiting each directed bond exactly once in a cycle. This 
is not the case on an infinite lattice, where the walker 
always finds new bonds which are not visited earlier and 
the number of visited sites keeps on growing. It was noted 
already that in 2d, the diameter of the region visited by 
the walker grows as V 1 / 3 , but the asymptotic shape was 
not investigated. 

The EW with noise is defined as follows: at each time 
step the walker rotates the arrow at its location by 0°, 
90°, 180° or 270° with probability e/3, 1 - e, e/3 or e/3 
respectively. We will show that the diameter of the region 
visited grows as N 1 / 2 for nonzero e and as N 1 / 3 when 
e = 0. 



III. NUMERICAL SIMULATIONS FOR 
EVOLUTION WITHOUT NOISE 

First we discuss the case e = 0. We denote the number 
of times different sites visited by the walker after TV steps 
by ri]y(x.) and the walker's square displacement from the 
origin by R 2 N . We evaluate tTn(x) (the over line rep- 
resents averaging over initial conditions), the variance 
of n(x) denoted by Var [tin(x.)], and the mean square 
displacement R 2 N by averaging over 10 6 different initial 
configurations. 



A. Mean square displacement 

According to the heuristic argument given in Ref. [||, 
if at time t the number of sites visited by the walker is 
S{t) , then in the previous 4S'(t) time steps, most of these 
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sites have been visited exactly 4 times except a small 
fraction at the boundary. As the cluster is seen to have 
few holes, it is nearly compact, and S(t) ~ D 2 (t), where 
D(t) is the diameter of the cluster, at time t. Thus we 
get 
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which implies that after TV steps, 
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Figure [3] shows the mean square displacement of the 
EW as a function of its length TV. The averaging is done 
over 10 6 realizations. The straight line, which is the best 
fit to the data has a slope 0.33 ± 0.01, consistent with 
Eq. ©. 



B. Average number of visits 

As seen in Fig. [TJ the cluster of visited sites is quite 
irregular in shape. Also sites that have been visited at 
least n times have rough boundaries with several islands 
of sites that have been visited fewer number of times than 
all the surrounding sites. However, if we average over 
different realizations of the initial arrow configuration, 
some interesting regularities are seen. 

In Fig. QJ we have plotted lines of njv(x) = for dif- 
ferent £ as indicated in the figure, for an EW of length 
TV = 10 6 averaged over 10 7 realizations. To obtain 
these lines, we add a diagonal bond between (x, y) and 
(x + l,y + 1) for each (x, y), and extend the definition 
of njv(x) to all real x by linear interpolation within each 
small triangle. The plot shows that these lines are nearly 
perfect circular in shape. 

The shape of rings, for large TV, can be defined by a 
function 



f(0) = lim 



TVi/3 : 



(3) 



where r/v(0) (0 < 9 < 2ir) is the angle dependent radius. 
If the shape is a perfect circle /(0) = constant, otherwise 
/(0) is a periodic function of that can be expressed in 
terms of Fourier cosine series 



f(6) = V a 4m cos(4m0). 



(4) 



Since the shape has fourfold symmetry the series will only 
have terms with m = Au (u = 0, 1, 2 . . . ). The vanishing 
of aim's for all m ^ then implies a circular shape. We 
define 
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J2j "jv(xj) cos(40j) 



(5) 
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FIG. 4: Contours of the constant average number of visits, 
rijv(x) for various values as indicated, for the EW of TV = 10 6 
steps averaged over 10 7 samples. 



as the normalized amplitude of the fourth Fourier mode. 
The summation j is over all the lattice points whose Eu- 
clidean distance from the origin lies between r and r + 1 . 
Here 6j is the angle that the vector Xj makes with the 
x-axis. This function, for a fixed r, has a well defined 
limit for N — » oo. In Fig.[5ja), we have shown A^(r) as a 
function of r for TV = 10 6 steps. The plot shows that for 
large r, A^(r) approaches zero with fairly large irregular- 
looking fluctuations. These fluctuations, for a fixed r, do 
not become smaller by statistical averaging or larger TV. 
These fluctuations occur because the lattice points lying 
between radii r and r + 1 are not distributed perfectly 
evenly along the ring. They are of number-theoretic ori- 
gin, and have been studied in the mathematics literature 
under the name of the 'Gauss circle problem' (l7| . The 
analysis of moments of A4 is therefore not very useful 
to estimate the shape and we have to adopt some other 
procedure. 

The simplest way to have a quantitative estimate of 
shape of rings is to fit the data and obtain the mean 
radius for various rings of constant njv( x )- In Fig- \Eft>)- 
we have shown such a ring and its fitting for njv(x) = 1 . 
The best fit gives the mean radius (R) = 145.436 ± 0.003 
for TV = 10 6 . For other rings also we find the error bars of 
the same order showing that the line of constant average 
number of visits are circular in shape within an error bar 
of 0.002%. The inset shows a close up of a particular 
region of the ring. 

We also calculate the root mean square deviation of 
distance, Ar(£), of points on the line of constant W^{x) = 
£ to the origin with mean radius {R)q- This is shown in 
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FIG. 5: (Color online) (a) A 4 (r) vs r for N = 10 6 steps, 
(b) Contour of rTjv(x) = 1 and its fit. The inset zooms a 
particular portion of the curve, (c) Ar(£)/ {R)( as a function 
of 9 for various rijv (x) = £ values. 



Fig. [5jc) as a function of 9 for various £. The plot shows 
that, Ar(C)/(i?)^ is of the order 10 -4 and decreases as 
{R)( increases. 

Another way to estimate the shape of the cluster 
formed by visited sites is to obtain various moments of 
the data. Since all the four directions are equivalent for 
the walker, we expect that njv(x) has a fourfold symme- 
try. For a given length N, we calculate (a; 4 ), (y 4 ) and 
(x 2 y 2 ) moments. If the shape of the cluster is perfect 



circular we would have 

(x 4 ) 



(y 4 ) 

(x 2 y 2 ) {x 2 y 2 ) 



(6) 



For N = 10 6 steps averaged over 10 6 initial realizations, 
we find that {x 4 )/(x 2 y 2 ) = (y 4 )/{x 2 y 2 ) = 3.007, which 
is consistent with the asymptotic value 3 deviation being 
only about 0.2%. 



C. Scaling of njy(x) 

For large N, rTjv(x), the average number of visits to 
the site x, depends only on |x|. Therefore, we expect 
that tTjv(x), satisfies the scaling form 



njv(|x|) = aN 1/3 F 



( |x| 



\bN 4 / 3 



(7) 



where F(y) is the scaling function. 

The scaling function can be determined as follows: Let 
xi and x 2 (|x 2 | > |xi|) be the distances of two different 
sites from the origin. The walker would have made sev- 
eral visits to xi before it first reaches x 2 . Afterwards, 
because of the local Eulcr-like organization of the arrows, 
both sites are visited equally often. Therefore, the differ- 
ence between the number of times sites at distances xi 
and X2 are visited remains bounded as N — > oo, i.e., 

nAr(xi) — 7T/v(x 2 ) = aN 1 / 3 [F(yi) — F(y 2 )] = constant. 

(8) 

This implies that F(y) must be a linear function of y. 
Using the freedom of choice of constants a, and 6, we can 
set F(0) = 1 and F{1) = 0. Therefore, we have 



F(y) = 



1-y for < y < 1, 
otherwise. 



(9) 



This simple form of F (y) was already noted by Toth and 
Veto for the problem in one dimension [l5|, [Ty| . The nor- 
malization condition J njf(x)dx = N gives ab 2 = 3/ir. 

In Fig. El we have plotted the finite size scaling of 
nW(xe x ) for N = 10 6 ,10 7 and 10 8 steps with a = 0.5 
and b = 1.38. In the same plot, we have also shown the 
scaling form given by Eq. ([9]). The plot shows that as 
N is increased, the scaled data approaches rather slowly 
towards the scaling form. In particular for smaller |x|, 
the approach to the asymptotic curve seems to be slow. 



D. Variance of njv(x) 

We also monitored the variance of n^ (x) as a function 
of N for different |x|. This is plotted in Fig. [7] The 
graph shows that Uar[njv(x)] increases slowly with N 
and suggests that yar[njv(x)] remains finite for all fixed 
x. This can be understood as follows: the variance of 
njv(x) arises only from the randomness in the initial visits 
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FIG. 6: (Color online)Data collapse of Wn(xb x ) for TV — 
10 6 , 10 7 and 10 s steps. The thick solid line is the scaling 
function. 

of walker to x. Once the local bonds have been organized 
into a near Euler circuit, the subsequent increments in 
njv(x) are nearly deterministic. 



E. Roughness exponent 

We define the surface of the set of visited sites as all 
visited sites that have at least one unvisited neighbor. 
Let W N denotes the variance of the distance from the 
origin of randomly picked surface formed by the EW of 
length TV. Then, W N is defined by 



^^ElN-W) 2 , (io) 

i=l 

where M is the number of surface points |xj| of the clus- 
ter, (R) is the average distance of a perimeter site, and 
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FIG. 8: (Color online) The width, Wn, of the boundary of 
the cluster formed by the EW vs its length TV. The data is 
averaged over 10 4 samples. The straight line is the best fit to 
the data with slope 8 = 0.136 ± 0.002. 

the overbar denotes averaging over different clusters. We 
define the width of the surface by square root of W N . 

In Fig. [51 we have plotted Wn for various TV, in a log- 
log scale, averaged over 10 4 different clusters. We observe 
that Wn ~ N s ~ R 3S ~ R a with 6 = 0.136 ± 0.02, 
which gives the roughness exponent a = 0.40 ± 0.06. 
The effective value of a seems to decrease with TV, and it 
is difficult to estimate its limiting value. It is consistent 
with the asymptotic value 1/3 expected for the Kardar- 
Parisi-Zhang (KPZ) surface growth process [HI]. Note 
that in order to think of growth of cluster of visited sites 
as a local growth process, we have to redefine time so 
that the radius of the cluster grows linearly in the new 
variable. 



IV. NUMERICAL SIMULATIONS FOR 
EVOLUTION WITH NOISE 
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FIG. 7: (Color online) Var [njv(x)] as a function of length TV 
for various |x|. The lines are the best fit to the data. 



For the evolution with noise also, we monitored R 2 N 
and njv(x) for various N by averaging over 10 6 different 
initial configurations. 



A. Mean square displacement 

For zero noise, for TV < 10 3 , R 2 N increases roughly 
linearly with N and then as TV 2 / 3 . In Fig. [31 we have 
also shown R? N as a function of TV for noise strengths 
e = 0.001, 0.01, and 0.1. For e = 0.001, there is only 
a small change from the Eulerian like behaviour. How- 
ever, when the noise strength is increased, there is a clear 
crossover seen in R 2 N from the Eulerian like to a simple 
random walk behaviour, i.e. R 2 N ~ TV. This crossover can 
be observed even for noise strength as small as (e = 0.01). 
That presence of noise changes the critical behaviour is 
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FIG. 9: (Color online) R 2 N {m)/R 2 N {0) vs m/N A with A 
0.74 ± 0.01 for various N. 




FIG. 10: (Color online) Contour of rw(x) = 1 for the EW 
of N — 10 5 steps in the presence of noise (e = 0.01). The 
averaging is done over 10 6 samples. The line is the best fit to 
the data. 



not so unexpected. In equilibrium critical phenomena, 
the well known Harris criterion (l9| characterizes a large 
class of systems where the critical behavior is substan- 
tially altered by the presence of disorder. In SOC mod- 
els, the Manna model [2(| with stochastic toppling rules 
is in a different universality class than the model with de- 
terministic toppling rules (e.g., the Bak-Tang-Wicsenfeld 
model) [H]]. I n fact, different types of stochasticity can 
yield different universality classes. For example, in the 
directed sandpile with stickiness one gets a different be- 
haviour than the stochastic Manna model (for the undi- 
rected case the situation is less clear [22|, HU ) . It is known 
that in some sandpile models, adding noise can change 
the transition from first order to continuous [24| . 

Let m — Ne be the number of 'm istak es' the walker 
makes in a walk of N steps and let R 2 N (m) represents 
the mean square displacement of a walk with m mis- 
takes. We would like to know how R 2 N {m) increases with 
to. Consider the effect of a single mistake. The rules 
of the Eulerian walk are such that as the walk evolves, 
the initial random arrow-directions are rearranged into 
an Euler circuit gradually. A wrong step would disrupt 
an evolving local Eulerian circuit. This local defect can 
be repaired, and a possibly different locally Euler-like cir- 
cuit formed, when the walker revisits the site. Therefore, 
there is not much change in R 2 N (m) for small to. How- 
ever, if to is larger (i.e., of the order TV), the walker keeps 
on making new mistakes before the old mistakes can be 
corrected, and the self-organization into Eulerian circuits 
is lost. This suggests that there should be a value to*, 
such that walks are nearly Eulerian for m < to*, and 
random walk-like for to > m* . We find tha t to* ~ 7V A 
with A < 1. Furthermore, we find that R 2 N (m) satisfies 
a scaling relation 

^(m)=^(0)G(jy , (II) 

where G(x) is the scaling function. In Fig. [5J we have 
plotted R%(m)/R%(0) vs m/N A for N = 1 x 10 5 , 5 x 10 5 



and 1 x 10 6 steps. A good collapse is obtained for 
A = 0.74 ±0.01. Equation (JTTJl can also be written 
as (using eN = m) ~R%{e) = ~R%(0)G (eiV 1 ^) . Then 
we have R%(0) ~ N 2/3 and R%(e) ~ TV, which implies 
that G{x) should increase as x 1 / 3 ^ 1- ^ and the cross over 
value N*(e) ~ e _1 / 3 ( 1-A ). A direct reliable estimate of 
N*(e) is not possible from our data and this scaling the- 
ory prediction is difficult to check from our simulations. 



B. Average number of visits 

In Fig. [TUl we have shown the contour of rTjv(x) = 1 
for an EW of length N = 10 6 steps in the presence of 
noise with strength e = 0.01. The averaging is done over 
10 6 initial conditions. The best fit to the data gives the 
mean radius (R) = 296.401 ± 0.015, i.e. circular within 
an error bar of 0.005%. The inset shows the close up of 
a particular region of the ring. On comparing this with 
the e = case, we see that the mean radius of an EW 
with 1% noise strength is about two times greater than 
the mean radius without noise. 

In Fig.rrjJ we have plotted rTN~(xe x ) vs x/N 1 ^ 2 , where 
e x is the unit vector along x-axis, for an EW in the pres- 
ence of noise of strength e = 0.01 for N = 1 x 10 5 , 1 x 10 6 
and 5 x 10 6 steps. The plot shows that, except near the 
origin where there is a log(iV) dependence, the data for 
various N collapse on top of each other. 

We also obtain the width, Wn, of the surface formed 
by the EW in the presence of noise. We find that 
Wn ~ N 1 / 2 ~ R as expected for random walks. Hence 
for nonzero e, the asymptotic shape of the cluster is not 
circular in individual realizations. Some circular symme- 
try is only seen in the ensemble averages (i.e., lines of 
constant n(x)). 
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FIG. 11: (Color online) Data collapse of riN(xe x ) vs x/N 1 ^ 2 
for the EW of N = 1 x 10 5 , 1 x 10 6 and 5 x 10 6 steps with 
noise strength e = 0.01. 

V. CONCLUSIONS 

We have studied an Eulcrian walker on a 2 dimensional 
square lattice using Monte Carlo simulations. In the ab- 
sence of noise, the mean square displacement (R%) ~ 



TV 2 / 3 . We find that lines of constant average number of 
visits seem to be perfectly circular. This result is not 
entirely unexpected, given the fact that Eulcrian walkers 
can be considered as a particular type of derandomized 
random walkers @. However, note that the result does 
depend on the fact that Eulerian walkers in two dimen- 
sion return to origin infinitely often, and while the av- 
eraged cluster shape would be expected to show circular 
symmetry, one does not expect each individual cluster to 
show spherical shape in dimensions d > 2. 

We also estimated the roughness exponent for the 
boundary of the visited region. This has a slow conver- 
gence to its asymptotic value, but our data is consistent 
with it belonging to the KPZ universality class. We also 
find that even a small randomness in the rule for walker's 
next step changes the asymptotic properties: the mean 
square displacement shows a crossover from the Eulerian 
like (R%) ~ N 2 / 3, to a simple random walk behaviour 
(R 2 N ) ~ N. In higher dimensions d > 2, an Eulcrian 
walker does not return to previously visited sites often, 
and one would expect a random walk like behavior even 
in the zero-noise case. 

DD thanks B. Toth and B. Veto for discussions. We 
thank J. Propp for a critical reading of the earlier version 
of the paper. 
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